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Abstract 

The multi-layer S = ~ square lattice Heisenberg antiferromagnet with up to 
6 layers is studied via various series expansions. For the systems with an odd 
number of coupled planes, the ground-state energy, staggered magnetization, 
and triplet excitation spectra are calculated via two different Ising expansions. 
The systems are found to have long range Neel order and gapless excitations 
for all ratios of interlayer to intralayer couplings, as for the single-layer system. 
For the systems with an even number of coupled planes, there is a second 
order transition point separating the gapless Neel phase and gapped quantum 
disordered spin liquid phase, and the critical points are located via expansions 
in the interlayer exchange coupling. This transition point is found to vary 
about inversely as the number of layers. The triplet excitation spectra are 
also computed, and at the critical point the normalized spectra appear to 
follow a universal function, independent of number of layers. 

PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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I. INTRODUCTION 



Low- dimensional quantum antiferromagnets exhibit many remarkable properties, and 
the study of these systems has been the subject of intense theoretical and experimental 
research in recent years. It is by now well established that one- dimensional Heisenberg 
antiferromagnets with integer spin have a gap in the excitation spectrum, whereas those 
with half-integer spin have gapless excitations. The former have a finite correlation length, 
while for the latter it is infinite with the spin-spin correlation function decaying to zero as 
a power law. In 2-dimensions, the unfrustrated square-lattice Heisenberg model has long 
range Neel order in the ground state. It has gapless Goldstone modes as expected. For 
the 3-dimensional Heisenberg model, one expects a stronger long-range Neel order |I] in the 
ground state. 

In recent years much interest has focussed on systems with intermediate dimensionality 
and on questions of crossovers between d = 1 and d = 2. The investigations @ showed that 
the crossover from the single Heisenberg chain to the two-dimensional antiferromagnetic 
square lattice, obtained by assembling chains to form a "ladder" of increasing width, is far 
from smooth. Heisenberg ladders with an even number of legs (chains) show a completely 
different behavior from odd-leg-ladders. While even-leg-ladders have a spin gap and short 
range correlations, odd-leg-ladders have no gap and power-law correlations. 

What happens to the crossover from the two-dimensional antiferromagnetic square lat- 
tice to the three-dimensional antiferromagnetic simple cubic lattice, obtained by assembling 
planes to form a multi-layer system? The spin-S" bilayer Heisenberg antiferromagnet has 
been well studied, and it has a second order transition |3|||] separating the Neel phase and 
dimer phase, and this critical points can be remarkably well fitted by the following linear 
relation: 



The systems with more than two planes have not been studied, as far as we are aware. As 
in the case with even and odd numbers of coupled chains, we may expect the Heisenberg 
antiferromagnet on an odd number of coupled planes to show a fundamentally different 
behavior from that for an even number of coupled planes. The multilayer Heisenberg anti- 
ferromagnet has recently attracted attention also due to its relevance to the understanding 
of the magnetic properties of cuprate superconductors (such as YBa 2 Cu 3 6+x .) containing 
two or more weakly coupled Cu0 2 layers |§-0. 

In this paper, we study the zero temperature properties of the n-layer, S — |, square lat- 
tice Heisenberg antiferromagnet, where each layer is composed of a nearest-neighbor Heisen- 
berg model on a square lattice, and there is a further antiferromagnetic coupling between 
corresponding sites of each layer. The system can be described by the following Hamiltonian: 
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S a j is a S=l/2 spin operator at site i of the ath layer, and denotes a pair of nearest 
neighbor sites on a square lattice. Here J\ is the interaction between nearest neighbor spins 
in one plane, and J 2 is the interaction between adjacent spins from different layers. We 
denote the ratio of couplings as y, that is, y — Jij J\. In the present paper, we study only 
the case of antiferromagnetic coupling, that is, both J\ and J 2 are positive. In the small 
J2/J1 limit, the model describes n weakly interacting 2D Heisenberg antiferromagnets, where 
each of them is Neel ordered and possesses gapless Goldstone excitations. While in the large 
J2/J1 limit, the system reduces to weakly coupled rungs, and each rung is described by the 
Hamiltonian 



The ground state of h rung is an S = singlet for even n, or an S = 1/2 doublet for odd n. 
This implies that there are some fundamental differences between even and odd numbers of 
coupled planes. For an odd number of layers, the system can be mapped to the well-known 
single layer Heisenberg antiferromagnet at the large J 2 limit, and so the system should have 
Neel order for all ratios of interlayer to intralayer couplings. For an even number of layers, 
on the other hand, the system has a definite gap at the large J 2 limit, so we expect there is 
a transition separating the gapless Neel phase and gapped quantum disordered spin liquid 
phase at a certain critical ratio of J2/J1, and this transition should lie in the universality 
class of the classical d = 3 Heisenberg model, as for the bilayer system M . 

To confirm the above arguments, we study these systems with up to 6 layers via various 
series expansions. This paper is organized as follows. In Section II, we develop two different 
Ising expansions for the 3 and 5-layer systems. The first one is the usual naive Ising expansion 
where we take the ^-component of the Hamiltonian as the unperturbed Hamiltonian, and 
the rest of it as the perturbation. In the second Ising expansion, we choose the eigenstates of 
^rung, which has an S = 1/2 doublet as the ground state, as the basis states. At the large J 2 
limit, H pi ane can be mapped to the well known spin-| single layer Heisenberg antiferromagnet 
described by an effective Hamiltonian H e ff, so we can take H Iung and the ^-component of 
if er r as the unperturbed Hamiltonian, and the rest as the perturbation. The first Ising 
expansion works best in the small J 2 / J\ region, while the second one works best in the large 
J 2 /Ji region. Series are calculated for the ground-state energy, staggered magnetization, 
and triplet excitation spectra. Extrapolating these series to the isotropic limit, we find 
the system has long range Neel order and gapless excitation for all ratios of interlayer to 
intralayer couplings, as for the square lattice. In Section III, the 4 and 6-layer systems are 
studied by means of expansion in the interlayer coupling J 2 , and the transition points are 
located. The last section is devoted to a summary. 

II. SYSTEM WITH ODD NUMBER OF COUPLED PLANES 

For n = 3, the eigenstates for ^ rung consist of 
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1) one Stot = 1/2 doublet with energy El nns = — 1 and eigenstates 

SL = ~\ ■ \{\, -\ -i)) = ^(1 lit) + I Til) - 2| iU» 



(5) 
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where the arrows represent, from left to right, the ^-components of S a (a = 1,2, ■■■,n), 
and the eigen states are denoted by (S tat , S* ot , E). S tot , S* ot and E are the total spin, the 
^-component of the total spin, and its corresponding eigenvalue, respectively. This is the 
ground state of h rung . 

2) another S tot = 1/2 doublet with energy E[ uns = and eigenstates 



.1 1 1 



% = l(o.--.°)> = -s(IUT>-|TU» 



(6) 



^t = |: 1(5,5,0)) = -L(| ITT) -ITU)) 

3) one S'tot = § quartet with energy E™ n& = 1/2 and eigenstates 



'-'tot — 2 ' Ko' 2'2^ _ ^ III) 
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SL = l- l(|^)) = ^(IITT) + ITIT)ITTD) 

3 3 1 

^tot = 2 : K2' 2' 2 = ^ 

So in the limit of J1/J2 = 0, the ground state of the whole system is the direct product 
of the ground states of each rung, which are degenerate with S = \ and S z = ±|, and so it 
is 2 N -io\d degenerate (N is no. of sites). A finite value of J\ lifts this degeneracy (to 2). In 
this 2 7V -dimensional subspace, one can define an effective Hamiltonian H e s which includes 
all interactions H p \^ ne . To first order in J%/ J 2 we obtain || 

-1 eff o off 



with effective coupling 



H e s — J e fi ^2 S° • (8) 



Jcs = Ji (9) 



Note that J e ff for Heisenberg ladder systems has been computed by N. Hatano and Y. 
Nishiyama ||, and their results can also be applied to the multi-plane Heisenberg antifer- 
romagnet discussed here. 
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Therefore the ground state energy per site to first order in J1/J2 will be 

E /N=(J 2 E^ + E?J eS )/n (10) 

where E^ 1 = —0.6693 is the ground state energy per site for the square lattice [[|. 

To find out the staggered magnetization at large J 2 limit, one needs to look at the results 
of the operator Sf — S 2 + Sf on the states |(|, — ~, —1)) and |(|, |, —1)): 

(St - SI + fig)|(i -\, -1)) = - -i -1)) + 0.942809|(|, -\, ±)> 

(11) 

(5? - 55 + ^DK^, ~, -1)) = i -1)) - 0.942809|(2 ~ i)> 

If we only keep the states which are the ground states of h Iung , and compare the above 
relations with that for the single layer square lattice, we can get the staggered magnetization 
at the large J 2 limit as 

1 5M sq 
M = ^E((-l) i+a ^)o = -oT- = 0-1706(6) (12) 

where M sq = 0.307(1) is the staggered magnetization for the single layer square lattice. || 
For other odd n values, the ground states of h rung are also two spin-| doublet states. 
In Table I, we listed the ground state energy E™ ng of h rnng , the effective coupling J e s and 
the staggered magnetization M at the large J 2 limit for n — 3,5, 7, 9, 11, and these are also 
shown in Figure 1, where we can see that as n increases, J e s/Ji is generally close to unity 
but increases monotonically, while M/M sq decreases about as rT 1 ! 2 . At the limit of n — > 00, 
we expect that J e s/Ji approaches a constant ||, while M/M sq decreases to zero. 

Since at both large and small y limits, the system can be reduced to a single plane 
Heisenberg antiferromagnet, we expect that unlike the bilayer Heisenberg antiferromagnets, 
there is no transition at any value of y. 

To study the properties of these systems, we develop following two different Ising expan- 
sions. 

The first one is the usual naive Ising expansion where we take the z-component of the 
Hamiltonian as the unperturbed Hamiltonian Hq, and the rest of it as the perturbation V: 
we write the Hamiltonian for the Heisenberg-Ising model as: 

H/J 1 = H + xV (13) 

where 

n n— 1 

a=l a=l i a,i 

(14) 

n n—l 

V = X] yi(Sn.jS%.-j + SajSaj) + y ^ y^XSnjSn^ A + ^a.i^a+lji) ~~ ^ e Qi,iSq,i 
a= l a=l i a,i 
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and £ a} i = ±1 on the two sublattices. The last term in both H and V is a local staggered field 
term, which can be included to improve convergence. The limits x = and x = 1 correspond 
to the Ising model and the isotropic Heisenberg model respectively. The operator H Q is taken 
as the unperturbed Hamiltonian, with the unperturbed ground state being the usual Neel 
state. The operator V is treated as a perturbation. It flips a pair of spins on neighbouring 
sites. The linked-cluster expansion method has been previously reviewed in several articles 
P^-fTS!, and will not be repeated here. 

The Ising series have been calculated for the ground state energy per site E /N, the 
staggered magnetization M, and the lowest lying triplet excitation spectrum A(k x ,k y ) for 
several ratios of couplings and (simultaneously) for several values of t up to order x 11 for 
n = 3 or order x 9 for n = 5. The triplet excitation spectrum is computed for n = 3 only. 
The resulting series for Eq/N and M for y = 0.5, 1, 2 are listed in Tables II, and the other 
series are available on request. 

There are three bands of the spin- wave dispersion for 3-layer systems. In the Ising limit, 
2 bands have initial excitations located in the side planes, and the third band has it in the 
middle plane. From the series one can see that all three bands have the following symmetry: 

A(k x , k y ) = A(k y , k x ) = A(n - k x , n - k y ) (15) 

Since at the large J 2 limit the system can be mapped to the well known spin-| single 
layer Heisenberg antiferromagnet with effective Hamiltonian H e g defined in Eq. (|8|), we can 
develop, in the basis of the eigenstates of h rung , another Ising expansion, the so called "mod- 
ified Ising expansion" , where we take H Tung and the ^-component of H efi as the unperturbed 
Hamiltonian, and the rest as the perturbation. That is, we write the Hamiltonian for the 
Heisenberg-Ising model as: 

H = H' + xV (16) 

where 

Tji J tt 1 7 rrcff,2 ncff,z . , V~~* otot,z 

U Q — J2-n r ung + ^eff / y Jj Oj + t CjOj 

(17) 

t/7 T tt T neS,z neS,z , rttat,z 

V — ^1-nplane — ^eff 2^ ~ 1 €i ^i 

5'j tot ' 2 is ^-component of the total spin for ith rung, and again the last term in both H' Q and V 
is a local staggered field term, which can be included to improve convergence. The limits x = 
and x = 1 correspond to the Ising model and the isotropic Heisenberg model respectively. 
The operator H' Q is taken as the unperturbed Hamiltonian, with the unperturbed ground 
state being the usual Neel state. The operator V is treated as a perturbation. The detailed 
eigenvalues and eigenstates of h rung and also the matrix elements V under the eigenstates 
of h Tnng for both n = 3 and n = 5 are available from the author. 

Here the Ising series have been calculated for the ground state energy per site Eq/N, 
and the staggered magnetization M for several ratios of couplings and (simultaneously) for 
several values of t up to order x 6 for n = 3 or order x 4 for n — 5. The resulting series for 
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Eq/N and M for y = 5, 10, 20 are listed in Tables II, and the other series are available on 
request. 

Having obtained the series for the above two Ising expansions we try to extrapolate the 
series to the isotropic point (x = 1). For this purpose, we first transform the series to a new 
variable 

5 = 1 - (1 - x) 1/2 (18) 

to remove the singularity at x = 1 predicted by the spin-wave theory. This was first proposed 
by Huse |13j and was also used in our earlier work on the square lattice case ||. We then 



use both integrated first-order inhomogeneous differential approximants flTJ] and naive Pade 
approximants to extrapolate the series to the isotropic point S = 1 (x = 1). The results 
for the ground state energy per site Eq/N, and the staggered magnetization M are shown 
in Figs. 2-3. The results for J2 = presented in the figures are taken from our earlier 
work on the single-layer square lattice 0, while the results for J1/J2 = are taken from 
the mapping of the system on to the single layer system. The asymptotic behaviour, Eq. 



(|10|), of the ground state energy is also shown in Fig. 2 as dotted line, and they match 
on to those Ising expansions very well. The results from two Ising expansion also match 
each other very well in certain regions of coupling. We note that M first increases for small 
J%j Ji, passes through a maximum at about J2/J1 — 1, and then decreases for larger values 
of J 2 /J\ ■ The reason that in the case of small J 2 / J\ , the interlayer coupling enhances the 
antiferromagnetic long-range order is that the system acquires a weak three dimensionality 
and quantum fluctuations are suppressed. The spin-wave theory for the bilayer Heisenberg 
antiferromagnet predicts that M should be proportional to (J2/ Ji) 1 ^ 2 in the small J2/J1 
region, we can expect this will happen for any multi-layer system, and so in Figs. 2-3 we 
show Eq and M as a function of (J2/ Ji) 1 ^ 2 /[J2J Ji) 1 ^ 2 + 1] to exhibit this behaviour. From 
Fig. 3, we can see that in the small J2/J1 region, the 5-layer system has slightly stronger 
Neel order than the 3-layer system, while in the large J%j J\ region, the 3-layer system has 
stronger Neel order. 

Figs. 4-6 show the three bands of triplet excitation spectra of the 3-layer system obtained 
from the naive Ising expansion. From these graphs, we can see that all of the dispersion 
relations have a minimum located at (k x , k y ) = (0, 0) (or (71", tt) by symmetry), where two of 
these three bands have a definite gap as long as J 2 7^ 0, while the third (the outer symmetry 
band) is consistent with a gapless spectrum for all values of J2. 



III. SYSTEM WITH EVEN NUMBER OF COUPLED PLANES 

For even n, the ground state of h rung is a S = singlet. For example, the eigenstates of 
ing for n = 4 consist of 
1) one S to t = singlet with energy i?Q Ung = (—3 — 2\/3)/4 and eigenstate 

Ko, 0, £r s )> = ci(i am + 1 mi» - c 2 (\ mt> + 1 uu» + <*a nit) + 1 mi)) (19) 



where 



3 + 2^ 12 + 7^ 9 + 5^ 

ci = — , c 2 = — , c 3 = — (20) 

6V26 + 15V3 6V26 + 15V3 6 J26 + 15 \/3 
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This is the ground state of H Tung . 

2) one S'tot = 1 triplet with energy E™^ = (— 1 — 2-\/2)/4 and eigenstates 



%t = "I = 
SL = : 



|(i, -l, £r s )> = c 4 (| HIT) - I Till)) + c 5 (| ITU) - I UU» 
|(i, o, K" 1S )) = c 4 (| IITT) - I TTU» + csd TITI) - I ITIT)) 
|(i, l, i?r g )) = cd(| ITTT) - I TTTO) + c 5 (| TTIT) - I TITT)) 



(21) 



where 



c 4 



l + v^ 



2^10 + 7^' 5 2^10 + 7^ 
3) another S'tot = 1 triplet with energy i?™ 1 ^ = —1/4 and eigenstates: 



3 + 2^2 



(22) 



^tot — — 1 : 



'S'tot — 1 : 



:i, -i, Er s )) = \{\ inn + 1 tiii) - 1 mi) 

;i,o,£r s )) = ^(IT!!T)-UTT!)) 
;i,i,Er s )) = ^(IITTT) + ITTTI)-|TITT)- 



itU» 



(23) 



TTIT)) 



4) another S'tot = singlet with energy E T 3 ung = (—3 + 2y / 3)/4 and eigenstate 

|(0, o, £r s )) = cad IITT) + I TTU)) - ci(| ITIT) + I TITO) - c 3 (| TUT) + I ITTO) (24) 

5) another S to t = 1 triplet with energy E™ ng = (—1 + 2v^2)/4 and eigenstates: 



'S'tot — — 1 : 



Cz 
D tot 



1 : 



;i, -i, £r s )) = c 5 (\ hid - 1 tiid) + c 4 (\ an) - 1 itid) 
;i, o, Er g )) = C5 (i mi) - 1 tiid) + c 4 (\ an) - 1 itid) 

;i, l, Er g )) = c 5 (| ITTT) - I TTTD) + c 4 (| TITT) - I TTIT)) 



(25) 



6) one S to t = 2 quintuplet with energy El ung = 3/4 and eigenstates 



CZ 

D tot — 


-2 : 


1(2, -2, Er s )) 


1 J' J' J' J' ^ 


cz _ 
D tot — 


-1 : 


i(2,-i^r g )) 


= \{\ HIT) + IIITD + I ITU) + 1 Till)) 


oz 

D tot — 


: 


K2,o,Er s )) = 


^=(1 IITT) + 1 ITIT) + 1 TUT) + IITTD + I 


CZ _ 

D tot — 


1 : 


K2,i,Er s )) = 


^(1 ITTT) + 1 TITT) + 1 TTIT) + 1 TTTD) 


cz 

'-'tot — 


2 : 


l(2,2,£r s )) = 


1 TTTT) 



So at the limit of J1/J2 = 0, the ground state of the whole system is the direct product 
of the ground state of each rung. Unlike the case for odd value of n, there is not degeneracy 
in this ground state, and the system has a definite excitation gap. 

The operator Sf — Sf + S| — Sf on the ground state |(0, 0, El nng )) gives 



(S? - S z 2 + S z 3 - S 4 z )|(0,0,^ ung )) = -1.45728|(1, 0,Er ng )) + 0.6036258|(1, 0, £4^)) (27) 
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Again unlike the case for odd n, the right hand side of the above equation does not contain 
the ground state |(0, 0, E r ung )), so in the large J 2 limit, the staggered magnetization M is 
zero, that is, the system is in a quantum disordered spin liquid state. So we expect there is 
a transition separating the gapless Neel phase and gapped quantum disordered spin liquid 
phase at a certain critical ratio of J2/J1, and we expect that this is true for any even number 
of layers. 

To locate the transition point, we can develop the expansions around the large interlayer 
coupling limit, that is, we can construct an expansion in l/yhy treating the operator H Tnng 
as the unperturbed Hamiltonian and operator H p \ ane as a perturbation. We have developed 
the expansion for the T = ground state energy per site Eq/N, the antiferromagnetic 
susceptibility Xi an d the lowest lying triplet excitation spectrum A(k x ,k y ) up to a certain 
order (see Table III) for n = 4 and 6. The resulting power series in 1/y for Eq/N and x are 
presented in Table III. The series for the minimum energy gap m = A (77, 7r) is also listed in 
Table III. The series for the excitations spectra A(k x , k y ) are available on request. 

To determine the critical point, we construct Dlog Pade approximants Jl4| to the x an d 
m series. Since the transition should lie in the universality class of the classical d = 3 
Heisenberg model (our analysis also supports this), we expect that the critical index for x 
and m should be approximately 1.40 and 0.71, respectively. The exponent-biased Dlog Pade 
approximants give the critical point (J x / J 2 ) c = 0.213(5) for n = 4, and (Ji/ J 2 ) c = 0.145(10) 
for n = 6. The critical points (</i/J 2 ) c as function of 1/n are shown in Fig. 7 (the results for 
n = 2 are taken from a previous calculation Q), where we can see they can be well fitted 
by {J\l J-i) c = 0.913/n — 0.250/n 2 . So (Ji/J 2 ) c should be proportional to 1/n at the large n 
limit. 

The triplet excitation spectra A(k x , k y ) at the critical point for n = 4, 6 are illustrated 
in Fig. 8(a): the results for n = 2 are also displayed here for ease of comparison. We can 
see, as in the case of the bilayer Heisenberg antiferromagnet [|J , the direct sum to the series 
is indeed consistent with the Pade approximants that one can construct. The spectra have 
their maximum value A(0,0) at (k x ,k y ) = (0,0). In Fig. 8(b), we plot the function for 
A(k x , k y )/A(0, 0) along high-symmetry cuts through the Brillouin zone, and we can see this 
function appears to be an universal function, independent of the number of layers n. We 
have no independent theory to explain it at present, but argue this may be due to the fact 
they belong to the same universality class. 



IV. CONCLUSIONS 

We have presented and analyzed various perturbation series expansions for multi-layer, 
S = h, square lattice Heisenberg model with up to 6 coupled planes. Like Heisenberg 
ladder system, the Heisenberg antiferromagnet on an odd number of coupled planes shows 
a fundamentally different behavior from that for an even number of coupled planes. 

For an odd number of layers, the system can be mapped to the well-known single layer 
Heisenberg antiferromagnet at the large J 2 limit, and so the system should have Neel order 
for all ratios of interlayer to intralayer couplings: that is, it is similar to the case of the square 
lattice. The staggered magnetization M at the large J 2 / Ji limit is computed up to 11- layers, 
and we found M should decrease to zero as n~ 1//2 in the large n limit. For the systems with 
3 and 5-layers, two Ising expansions have been developed, and series are calculated for 
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the ground-state energy, the staggered magnetization, and the triplet excitation spectra. 
Extrapolating these series to the isotropic limit, we find the system has long range Neel 
order and gapless excitations for all ratios of interlayer to intralayer couplings. 

For an even numbers of layers, expansions in the interlayer coupling have been developed 
for both 4 and 6 layers. The systems turn out to have a second order transition separating 
a gapless Neel phase and a gapped quantum disordered spin liquid phase, where the critical 
point is determined to be {J\/J2)c = 0.213(5) for n — 4, and {J\/J2)c — 0.145(10) for 
n — 6, and this critical point should be proportional to 1/n it the large n limit. The 
triplet excitation spectra are also computed, and at the critical point, it turns out that the 
excitation spectrum, after normalization by its maximum value, appears to be an universal 
function, independent of the number of layers. This may be due to the fact that they belong 
to the same universality class. 
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FIGURES 



FIG. 1. The effective coupling J e g and the staggered Magnetization M for odd n at the 
limit J2/J1 — * 00 as function of n -1 / 2 , and the solid line is the least square polynomial fit 
M/M sq = 1.125n- 1 / 2 - 0.532n -1 + 0.431n" 3 / 2 . 

FIG. 2. The rescaled ground-state energy per site Eq/N for 3 and 5-layers as function of 
{J2/J1 ) 1 ^ 2 / [{J2/J1 ) 1 ^ 2 + !]• The solid (open) points with error bars are the estimates from the 
naive (modified) Ising expansions, and the dotted curve at large (J2/ Ji) 1 ^ 2 region is the asymptotic 
behaviour, i.e., Eq. JlC|). 

FIG. 3. The staggered magnetization M for 3 and 5-layers versus (J2/ Ji) 1 ^ 2 /[(J2/ Jl) 1 ' 2 + !]■ 
The solid (open) points with error bars are the estimates from the naive (modified) Ising expansions. 

FIG. 4. Plot of symmetric spin-triplet excitation (in side planes) spectrum A(k x ,k y ) (derived 
from the naive Ising expansions) along high-symmetry cuts through the Brillouin zone for the 
three-layer system with coupling ratios J2/J1 = 0,0.5,1 (shown in the figure from the bottom to 
the top respectively). 

FIG. 5. Plot of antisymmetric spin-triplet excitation (in side planes) spectrum A(k x ,k y ) (de- 
rived from the naive Ising expansions) along high-symmetry cuts through the Brillouin zone for the 
three-layer system with coupling ratios J2/J1 = 0, 0.5, 1, 2, 4 (shown in the figure from the bottom 
to the top respectively). 

FIG. 6. Plot of spin-triplet excitation (in middle plane) spectrum A(k x ,k y ) (derived from the 
naive Ising expansions) along high-symmetry cuts through the Brillouin zone for the three-layer 
system with coupling ratios J%j J\ = 0, 0.1, 0.25, 0.5, 1 (shown in the figure from the bottom to the 
top respectively). 

FIG. 7. Plot of the critical points (Ji/ ' J%)c versus 1/n for even number of coupled planes. The 
solid points with error bars are the estimates from expansions about interlayer coupling, and the 
solid line is the least square polynomial fit (J\/J2)c = 0.913/n — 0.250/n 2 . 

FIG. 8. Plot of the triplet excitation spectrum A(k x ,k y ) (a) and spectrum after normalization 
by its maximum value (b) along high-symmetry cuts through the Brillouin zone for the n = 2, 4, 6 
systems at the transition point, the lines are the estimates by direct sum to the dimer series, and 
the points are the estimates of the Pade approximants to the series. 
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TABLES 



TABLE I. the ground state energy E r nn& of h rung , the effective coupling J e g and the staggered 
magnetization M in the limit J2/J1 — ► 00 for odd n up to n = 11. 



n £ ( T s /n J cff /Ji M/M Bq 

3 -0.33333333 1 0.555556 

5 -0.38557725 1.016938829 0.435575 

7 -0.40517710 1.034391333 0.372626 

9 -0.41514686 1.050406957 0.331978 

11 -0.42109939 1.064819159 0.302836 



TABLE II. Series coefficients of naive and modified Ising expansions of three layer system for 
the ground-state energy per site Eq/{NJ\) and the staggered magnetization M, for J2/J1 = 0.5, 1, 2 
(for naive Ising expansion) 5, 10, 20 (for modified Ising expansion) and t = 0. Nonzero coefficients 
x % up to order i = 10 for naive Ising expansion or to order i = 6 for modified Ising expansion are 
listed. 



i 






E /(JiN) 










M 


















Naive Ising expansion 














J2/J1 = 0.5 


-1 


J2/J1 = 1 


1 


J2/J1 = 2 


J 2 /J! = 0.5 


-1 


J2/J1 = 1 


-1 


J2/J1 = 2 







-3.33333333x10" 


-4.16666667x10- 


-5.83333333X10- 1 


5.00000000x10- 


5.00000000x10- 


5.00000000x10- 


1 


2 


-1.46708683x10- 


-1 


-1.53703704x10- 


1 


-2.23809524X10" 1 


-7.98687122x10- 


-2 


-7.14609053x10- 


2 


-8.68027211x10- 


2 


4 


1.10044693x10- 


-3 


5.08576750x10- 


4 


1.36582616xl0~ 3 


-6.95556119x10- 


-3 


-5.95673893x10- 


3 


-8.37837732x10- 


3 


6 


-1.39063516x10- 


-3 


-1.97000255x10- 


3 


-2.95007447xl0" 3 


-4.80430415x10- 


-3 


-5.53318319x10- 


3 


-7.46816918x10- 


3 


8 


-5.11831870x10- 


-4 


-6.87771224x10- 


4 


-1.22506627xl0 -3 


-3.03361100x10- 


-3 


-3.35287068x10- 


3 


-5.27035727x10- 


3 


10 


-2.75677029x10- 


-4 


-4.09491760x10- 


4 


-6.02299585X10" 4 


-2.01622805x10- 


-3 


-2.42979325x10- 


3 


-3.42232873x10- 


3 












Modified Ising expansion 














J2/J1 = 5 




Ja/Ji = 10 




J2/J1 = 20 


J2/J1 = 5 




J 2 /Ji = 10 




J2/J1 = 20 







-1.83333333 




-3.50000000 




-6.83333333 


2.77777778x10- 


-1 


2.77777778x10- 


-1 


2.77777778x10- 


1 


1 


0.00000000 




0.00000000 




0.00000000 


6.97167756x10- 


-2 


3.70370370x10- 


-2 


1.91158901x10- 


2 


2 


-1.21179736x10- 


-1 


-9.07296550x10- 


2 


-7.38090200xl0~ 2 


-4.36538214x10" 


-2 


-5.67336122x10- 


2 


-6.04157972x10- 


2 


3 


-1.41804384x10- 


-2 


-3.48154003x10- 


3 


-6.42216094X10" 4 


-1.35350826x10- 


-2 


-7.15990062x10- 


3 


-3.64767685x10- 


3 


4 


4.43598485x10- 


-3 


3.25869698x10- 


3 


1.90188695xl0~ 3 


-5.47105059x10- 


-3 


-4.79712404x10- 


3 


-6.78984443x10- 


3 


5 


5.30206875x10- 


-3 


1.33866469x10- 


3 


3.36884895X10- 4 


9.81694291x10- 


-3 


2.83116612x10- 


-3 


5.48197309x10- 


4 


6 


-1.30627826x10- 


-3 


-1.28650672x10- 


3 


-8.82824083X10- 4 


-3.28072052x10- 


-3 


-5.39705234x10- 


3 


-5.29473180x10- 


3 



TABLE III. Series coefficients for expansions about interlayer coupling of the ground-state 



energy per site Eq/{NJ2), the (tt, it) gap m/ 'J2, and the antiferromagnet susceptibility \ f° r 4 and 
6 layer systems. Nonzero coefficients (Ji/^)* up to maximum order carried out are listed. 









4 layers 








6 layers 




i 


E /(NJ 2 ) 




m/J2 


X 


E /(NJ 2 ) 




m/J2 


X 





-4.040063509x10- 


1 


6.589186226xl0 _1 


1.699358737 


-4.155961890x10- 


1 


4.915817770xl0 _i 


2.403792997 


1 


0.000000000 




-2.149829914 


1.161082473X10 1 


0.000000000 




-2.222865048 


2.332835586 xlO 1 


2 


-4.415865088x10- 


1 


-9.251531400X10" 1 


6.622899086 xlO 1 


-4.924510549x10- 


1 


-1.864039733 


1.902981798 xlO 2 


3 


-2.554087280x10- 


1 


-6.521784931 


3.610950193X10 2 


-3.192744899x10- 


1 


-1.517689040X10 1 


1.492331708 xlO 3 


4 


-5.884519404x10- 


1 


-8.897766802 


1.876359629 x10 s 


-1.446337784 




-2.769539610X10 1 




5 


-9.377312694x10- 


2 


-2.575353400X10 1 


9.557326563 xlO 3 


-7.523886385x10- 


1 


-1.455619317xl0 2 




6 


-1.120389803 




-4.663997073 xlO 1 












7 


4.729204915x10- 


1 


-3.236835733 xlO 2 
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0.6 





(0,0) (tt,0) (tt.tt) (0,0) 



< 




(0,0) (tt,0) (tt.tt) (0,0) 



